Light transport in cold atoms and thermal decoherence 
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By using the coherent backscattering interference effect, we investigate experimentally and the- 
oretically how coherent transport of light inside a cold atomic vapour is affected by the residual 
motion of atomic scatterers. As the temperature of the atomic cloud increases, the interference 
contrast dramatically decreases emphazising the role of motion-induced decoherence for resonant 
scatterers even in the sub-Doppler regime of temperature. We derive analytical expressions for the 
corresponding coherence time. 
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Originally motivated by astrophysical purposes, wave 
transport in opaque media was first analyzed through a 
detailed balance of energy transfers within the scattering 
medium, leading to radiative transfer equations . One 
main physical ingredient of this theory is a random-phase 
assumption, thereby discarding all possible interference 
effects. In its simplest form, this theory predicts a diffu- 
sive transport with spatial diffusion constant D = i'^ /3t, 
where i is the scattering mean free path and r the trans- 
port time. However, it is now clear that interference can 
survive disorder average in the elastic regime and dra- 
matically alter wave transport, leading for example to 
the physics of weak and strong localization Q . There has 
thus been much effort in solid state physics to study, mas- 
ter and circumvent any possible phase-breaking mecha- 
nism to fully access the regime of coherent transport Q ■ 
Coherence loss phenomena can be characterized by a 
phase-breaking time r* and a corresponding coherence 
length L$ = y/Dr^ |2| beyond which interference ef- 
fects are essentially washed out. The coherent transport 
regime (also known as the mesoscopic regime) is then 
reached when the medium size L is smaller than For 
instance, in metals or semi-conductors, it prevails up to 
L$ sa 1 — 10/im even at very low temperatures and very 
pure samples, thus requiring the use of devices with mi- 
crometer size 0. 

During the last two decades, the field got a renewed at- 
tention in systems using electromagnetic waves and, 
among the large variety of scattering media employed so 
far, cold atomic vapors have recently emerged. An im- 
portant feature of these media is the sharp resonance of 
the scattering cross-section (of order on resonance, A 
being the wavelength), allowing for a continuous tuning 
of the scattering mean free path i via the light frequency. 
Phase-breaking mechanisms in cold atomic vapors can be 
efficiently probed by using the coherent backscattering ef- 
fect (CBS), a paradigmatic two- wave interference effect 
in multiple scattering . In the limit of low light inten- 
sity, an important mechanism has been identified 0|. It 



is rooted in the Zeeman degeneracies of the atomic in- 
ternal structure but could be healed by applying a mag- 
netic field to the sample (2|]. Another important phase- 
breaking mechanism is the residual thermal motion of 
atomic scatterers inside the cold vapor and is the main 
topic of this paper. Indeed, because of the Doppler and 
recoil effects, light frequency is changed upon scatter- 
ing and transport is no more elastic. In previous pa- 
pers |M 1 have shown how the incoherent transport 
is affected by such effects. For typical alkali magneto- 
optical traps (MOT), Doppler- induced frequency shift is 
the dominant effect. Not surprisingly, an important pa- 
rameter is the ratio of the typical Doppler shift kv (v 
being the ID rms velocity, and k — 2'k/X the incom- 
ing light wavevector) to the width F of the atomic reso- 
nance. When few 3> F, a single scattering event is enough 
to bring the photon completely out of resonance thereby 
altering transport significantly [inl |. When few <gC F, a 
single scattering event only slightly modifies the photon 
frequency and leads, as a cumulative effect, to diffusion 
in frequency space |9j- 

In this Letter, we report the first unambiguous exper- 
imental evidence of a phase-breaking mechanism in light 
transport induced by the residual atomic motion . We 
first describe the experimental procedure and our main 
result. We then give an analysis of the physical ingredi- 
ents at the heart of the light coherence reduction. We 
finally compare our data to the Monte-Carlo simulations 
described in |l2l| . appropriately modified to include the 
effect of the atomic velocity distribution ■ 

Our CBS experimental setup has been described in de- 
tails elsewhere the main difference here being the 
heating procedure to increase v. It is achieved by switch- 
ing off the magnetic gradient of the MOT and then ex- 
posing the cloud to a 200^s long, slightly red-detuned, 
optical molasses. An increasing amount of heating is 
obtained by tuning the molasse frequency closer to reso- 
nance. The resulting 3D velocity distribution is measured 
by releasing atoms from the trap and imaging their ballis- 
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FIG. 1: CBS enhancement factor (full circles) as a function 
of the ID typical atomic velocity v and its comparison to a 
Monte-Carlo calculation assuming either a Gaussian (dashed 
line) or a Lorentzian (solid line) atomic velocity distribution. 
The CBS peak is recorded in the helicity non-preserving chan- 
nel, with resonant laser light (5 = 0). The optical thickness 
of the cloud is set constant at 6 ~ 13. For Rubidium, the 
velocity scale T/k is 4.6 m/s. 



tic expansion. With some word of caution (one side-effect 
of this heating technique is to change the atomic velocity 
distribution which can differ from a Gaussian one), the 
temperature of the cloud is fc^T = mv'^. The number of 
atoms in the cloud can be adjusted independently of T so 
as to maintain a fixed optical thickness b — L/£. In Fig.Q] 
we observe a fast decrease of the CBS enhancement fac- 
tor as V is increased, the value obtained as w ^ being 
set by the Rubidium Zeeman degeneracy. In particular, 
the decay is faster than expected from the naive criterion 
kv M T. The measurements are performed in the helicity 
non-preserving channel with a resonant laser (detuning 
6 = Lu — ujQ = 0, where lo and luq are the laser and atomic 
angular frequencies). The measured angular width of the 
CBS peaks ~ l/ki {k£ ~ 1000 typically) varied by 
less than 10% in the course of the experiment, which 
proves that the reported decrease is not driven by the 
angular resolution of the apparatus. Before commenting 
Fig. n let us discuss the physical ingredients of motion- 
induced CBS reduction for the simple double scattering 
configuration shown in Fig. We assume here that the 
atomic motion is classical ^IJ]. Atoms, distributed in 
space with number density p <C A:^, are illuminated by 
a plane wave of angular frequency uj and wave vector k 
and build up a dilute effective medium characterized by 
its complex index of refraction n w 1 + pa/2, convolved 
by the atomic velocity distribution. The atomic com- 
plex polarizability is a = -SnTk-^ / {S + ir/ 2). The CBS 
signal can be understood as an interference between a 
"direct" scattering path (solid arrows) and its "reverse" 
counterpart (dashed arrows). To each path is associated 




FIG. 2: (color online) CBS dephasing induced by the atomic 
motion, in the simplest case of two atoms moving with dif- 
ferent velocities vi and V2. The CBS effect builds on the 
two-wave interference between amplitudes associated to re- 
versed multiple scattering paths (arrows) . The Doppler effect 
induces frequency redistribution of the scattered light. Al- 
though the outgoing frequencies are identical, the intermedi- 
ate frequencies differ, yielding a phase coherence loss and thus 
a reduction of the CBS contrast. 

a complex amplitude both incorporating scattering and 
propagation in the effective medium. For scatterers at 
rest, with no Zeeman degeneracies, these amplitudes are 
exactly balanced at backscattering (outgoing wave vector 
— k) yielding full CBS contrast. For Zeeman degenerate 
atoms, like Rubidium, we recover the results of :6|. Let 
us now consider the impact of the atomic motion. Be- 
cause of the Doppler effect, the light frequency in the 
laboratory reference frame is modified. For the direct 
path, it is = cij + (k' — k).vi where k' is the intermedi- 
ate wavevector. After scattering by atom 2, the outgoing 
frequency is w' = uj + (k' — k).vi — (k + k').V2. For 
the reverse path, the frequency between the two atoms 
is LOr = UJ — (k' -|- k).V2 7^ LOd, while the final frequency 
is again lo' . Thus, although the outgoing frequencies are 
identical and do interfere, the intermediate frequencies 
differ by an amount ~ kv, resulting in different complex 
amplitudes and reducing on average the CBS interference 
contrast. Similarly, the frequencies seen by the atoms in 
their rest frames are different along the direct and reverse 
paths, leading to different scattering phaseshifts. More 
specifically, for small velocities kv one can estimate 
the typical phase difference A$ induced by the atomic 
motion between the reversed paths: 

A$ = A$, + A$p ~ [d^<i>, + ked^n] kv « kv/T (1) 

The first term corresponds to the scattering phaseshift, 
while the second describes the phaseshift associated to 
propagation over a distance £ in the effective medium. 
Although approximate, eq. 1^ provides useful insights 
about the impact of S at finite v (see below) . For larger 
scattering orders {N > 2), the elementary scattering and 
propagation processes have to be chained and randomly 
distributed phase differences will accumulate, decreasing 
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further the CBS contrast, see FigO In order to accu- 
rately compute the shape and height of the CBS peak, 
one must average over all atomic degrees of freedom, both 
external (position and velocity) and internal (magnetic 
quantum number), but one must also incorporate the 
experimental geometry such as the shape, size, optical 
thickness of the cloud as well as the waist and linewidth 
of the CBS probe laser beam. A Monte-Carlo approach 
as described in is then necessary, where we now in- 
clude the effect of the atomic velocity distribution % As 
mentioned earlier, some of the measured velocity distri- 
butions exhibited clear deviations from a Gaussian with 
an excess of fast atoms. Such atoms are responsible for 
large Doppler effects and are thus likely to kill CBS more 
efficiently. We have used in our simulations either a 
purely Gaussian velocity distribution or a Lorentz-like 
distribution (where we define v — vq/ ('\/8Zn2), vq being 
the FWHM- velocity) . The results, displayed in Fig. ^ 
show a good qualitative agreement, with a fast decay 
of the enhancement factor with increasing v. However, 
the Gaussian distribution show deviations at low v while 
the Lorentz distribution gives a better agreement. This 
pinpoints the important role played by a small fraction 
of atoms with large velocities. A quantitative agree- 
ment would require a better knowledge and control of 
the initial atomic velocity distribution. As can be seen, 
the CBS contrast decreases faster than naively expected. 
This is because higher scattering orders contribute signif- 
icantly to the CBS effect. Indeed, at large optical thick- 
ness {b > 10), they contribute altogether to « 30% of 
the CBS peak height. As increasing scattering orders 
are more and more sensitive to motion-induced phase- 
breaking, this explains the reported observations. 

Since the CBS scrambling originates from the 
frequency-dependent response of the sample (which is 
maximum around resonance), one could expect to re- 
cover a full interference contrast by detuning the laser 
at |i5| » kv. This is an important issue for the prospec- 
tive CBS observation using a room-temperature vapor. 
Fig. El displays Monte-Carlo simulations of the CBS con- 
trast as a function of the rms velocity v at S = and 
S ^ kv. Calculations are made for N = 2 (squares) and 
iV = 40 (circles) in a semi-infinite medium. As can be 
seen (curves a and h, N = 40), at low velocities fci) <C F 
the interference contrast is nearly frequency independent. 
This behavior is due to a compensation between the vari- 
ations of the scattering and propagation terms. Indeed, 
the phase difference in eq. ^ is essentially proportional 
to the transport time r (bracketed term), which was 
shown to be frequency-independent near the resonance 
and equal to F"^ As mentioned earlier and expressed 
in eq. 10), the decay rate is faster for higher scattering 
orders. At larger velocity (curves 1 and 2, N — 2), 
Fig. El shows an increased frequency dependence of the 
interference contrast, which we interpret as an unperfect 
balance between the scattering and propagation phase- 
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FIG. 3: (color online) Numerical calculation of the CBS inter- 
ference decay for multiple scattering order N — 2 (squares) 
and = 40 (circles). The upper velocity range refers to 
curves a and b, while the lower one refers to curves 1 and 
2. Note the very different velocity ranges. We compare two 
situations: resonant excitation S — {1 and a) and detuned 
excitation 5 > fct; (5 = SOP for 2 and S = 5T for b). Al- 
though the CBS contrast is improved, a detuned excitation is 
not sufficient to circumvent the motion-induced decoherence. 



shifts in eq. The interference contrast is partially 

restored in the \S\ 3> kv limit, where the motion- induced 
decoherence is now essentially driven by the propagation 
term in eq. ^ of order kv/T. The interference thus still 
decreases towards zero with increasing v. 

Beyond the specific case of CBS, one can quantitatively 
characterize the motion-induced phase breaking by defin- 
ing the phase-breaking time t$ and its associated coher- 
ence length L$. We discuss here the simplest low- velocity 
case, kv ^ F. At each scattering event, the photon de- 
tuning is slightly modified and this cumulative process 
turns into a random walk at sufficiently large scattering 
orders 1. Its dispersion is « kv^/N which we also 

assume to be ^ F. If the atomic velocities are Gaussian 
distributed and independent, the averaging over the de- 
grees of freedom can be performed analytically, and the 
resulting interference contrast decreases like: 



cm{v) = exp 



/kv^ ^ 

u [y 



(2) 



Note that this quantity is independent of the initial de- 
tuning S as discussed before. This expression is also valid 
for larger velocities, in the off-resonant limit S 3> VN kv 
where the convolution by the velocity distribution can 
be neglected. The decay of the interference contrast de- 
scribed by eq. (j^J is not exponential with N, in contrast 
to other sources of decoherence like the atomic internal 
structure [^. Nevertheless, one can still define a criti- 
cal scattering order = 3(3fcu/2F)~^/^. This means 
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that the motion-induced phase-breaking mechanism can 
be characterized by: 
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Therefore, the mesoscopic condition L$ > L reads: 
kv 1 



(3a) 
(3b) 

(4) 



For our lowest temperature of 40 /iK {v — 6.5 cm/s), 
eq. (|3bp yields L$ ~ 4£ and condition Q is fulfilled for 
optical thicknesses b < A. 

At this point, an interesting comparison with electronic 
transport in metals can be made. There, the temper- 
ature dependence of r$ is governed by several differ- 
ent mechanisms: electron-electron scattering (oc T~^/^), 
electron-phonon scattering (oc T^^) and spin-flip scat- 
terin g (w eakly T-dependent) which dominates at low 
T UJS- In our case, t$ is set by the internal struc- 
ture [6j (independent of T) and by the motion-induced 
decoherence, which, according to eq. H3al) . yields t$ oc 
T-i/3_ At T = 40 fiK, T$ is mostly determined by the 
internal structure (see Fig. 

A simple dynamical picture of the motion-induced co- 
herence loss has been given by Golubentsev In this 
picture, phase-coherence is lost when the length of a mul- 
tiple scattering path varies (due to the motion of the 
scatterers) by A during the time taken by the light to 
follow this path. Golubentsev found expressions equiv- 
alent to ours, with however a huge difference: for the 
non-resonant case he considered, the transport time is 
T = £/c < 1 ps while for our narrow resonance this time 
is « 27 ns, more than 4 orders of magnitude higher. 
It is only because photons are "slowed down" inside the 
resonant medium that we are able to observe the destruc- 
tion of phase coherence. 

In conclusion, we reported the first observation of a 
motion-induced phase breaking of the CBS effect. The 
observed CBS reduction is due to the fast-varying fre- 
quency response of the medium in the vicinity of a res- 
onance, which enhances the impact of frequency redis- 
tribution caused by Doppler effect. The rapid decrease 
of the interference contrast indicates that the phase co- 
herence length L$ associated with the residual motion of 
atoms in our coldest sample is not very large compared 
to the scattering mean free path, and smaller than the 
actual size of the cloud. According to the mesoscopic cri- 
terion derived in this paper, extremely low temperatures. 



in the nK range, are probably required to fully access in- 
terference effects in optically-thick cold atomic gases. In 
this regime however, the recoil effect is expected to be- 
come the main source of decoherence. 
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